Dimensional collapse and fractal attractors of a system with fluctuating delay times 
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A frequently encountered situation in the study of delay systems is that the length of the delay 
time changes with time, which is of relevance in many fields such as optics, mechanical machining, 
biology or physiology. A characteristic feature of such systems is that the dimension of the system 
dynamics collapses due to the fluctuations of delay times. In consequence, the support of the long- 
trajectory attractors of this kind of systems is found being fractal in contrast to the fuzzy attractors 
in most random systems. 
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Introduction. In ubiquitous natural and laboratory sit- 
uations the action of time delayed signals is an essential 
ingredient to understand the system dynamics. For in- 
stance, in optical and physiological systems [l|43|, a finite 
transmission speed usually leads to delayed reactions to 
signals from remote origins and/or time shifts among sig- 
nals taking different paths. In mechanical engineering the 
rotation of the workpiece is the origin of time delay ef- 
fects and in feedback control methods time delayed 
signals are utilized to stabilize the system behavior [5] . 

A frequently encountered situation which was rarely 
addressed previously [6], is that the length of the de- 
lay time often varies with time. A biological example 
is that the reproduction cycle of animals fluctuates with 
the environment change. In mechanical engineering, vi- 
brations of the machine tool and/or the workpiece may 
also change the length of the delay time [8|. As a mod- 
ern control scheme, the time variation of delays has been 
found to enhance the control efficiency [9]. 

The relevance of fluctuating delay times for all these 
phenomena motivates us to gain a better understanding 
of its influence on the system dynamics. In this contribu- 
tion, we will start from a very simple map system where 
the delay time takes only the value of one or zero dis- 
crete time steps. High dimensional systems with longer 
fluctuating delays show similar behavior [111, [T2j. 

A characteristic feature of such systems is that the sys- 
tem dynamics collapses to a low dimensional subspace 
due to the fluctuation of delay times. One would expect 
that such a dimensional collapse will have some inter- 
esting consequences on the system dynamics. As an ex- 
ample, we show in this paper its influence on the fractal 
properties of the system attractors. 

Models. For simplicity, we focus on the map system 



ax. 



(1) 



where k and b are constants and the time dependent ran- 
dom variable ct takes only value or 1. The probability 
of ct taking the value 1 is denoted as p. For a parameter 
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setting with 6 = the length of the delay time of our sys- 
tem can vary randomly between one or zero time steps 
with the system dynamics simply switching between the 
Henon map (e^ = 1) and the logistic map (e^ = 0). To 
illustrate the importance of the delay time variation and 
the resulting dimensional collapse a comparison is made 
with the well-studied case of the random Henon map in 
[13] with b and ct being an uniform noise in the in- 
terval [0, 1], where the delay time is constantly 2. In the 
following the parameter a is fixed as a = 1.1 if it is not 
stated otherwise. 

Numerical results. As depicted in [13] to study the ef- 
fect of randomness on attractors of a system like Eq.([T]) 
one may consider two different situations, either the 
snapshot attractor formed by an ensemble of points at 
a single instant of time with all points starting from dif- 
ferent initial conditions and evolving under a same given 
realization of the external noise, or the long- trajectory at- 
tractor formed by a long trajectory segment starting from 
a single initial condition and evolving under a given re- 
alization of the external noise. The main issue addressed 
here will be the difference between the long-trajectory at- 
tractors of the random delay case and the random Henon 
case. 

It was shown [l3| that for the random Henon case 
6/0 the long-trajectory attractor is a fuzzy attractor 
with a smooth density of points while the correspond- 
ing snapshot attractor may be fractal. An example with 
b = 0.3, k = 0.05 is presented in Fig(T] (a) and (b). In 
contrast, for the random delay case b = we find that 
the long-trajectory attractor turns out to be fractal in 
consequence of the delay time variation and the dimen- 
sional collapse. In Fig(T] (c) and (d) the long-trajectory 
and snapshot attractors are shown for a random delay 
case with 6 = 0, /c = 0.3 and p = 0.5. Moreover, as 
demonstrated the long-trajectory attractors obtained are 
identical for two runs with different initial conditions and 
different realizations of the noise ct. This is different from 
previously investigated random systems where in general 
different realizations of external noise lead to different at- 
tracting sets lying in different regions of the phase space. 
The snapshot attractor of the random delay case is time 
dependent and turns out to be at probability 1 a one- 
dimensional curve (see Fig(T](d)). 



2 



To show quantitatively the difference in the fractal 
properties between the random delay case and the ran- 
dom Henon case we calculated the capacity and infor- 
mation dimension of the attractor via the box-counting 
algorithm. The variation of the number of non-empty 
boxes N{1) versus the box size / is plotted in FiglS By 
definition the slope of the log- log plot of N{1) versus l/l 
gives the capacity dimension Dq of the attractor. Four 
curves in Figl2] correspond to the four attractors shown 
in FigdJ To guide the eyes, three lines of slope 1, 2 and 
1.26 (the capacity dimension of the Henon attractor) are 
plotted in the same figure. One can easily see from the 
plot that the capacity dimension of the long-trajectory 
attractor of the random Henon case and the random de- 
lay case is very close to 2 and 1.26 respectively. This 
confirms our above observation from Fig(T]that the long 
trajectory attractor of the random Henon case is a fuzzy 
attractor with a smooth density of points while the long 
trajectory attractor of the random delay case is a fractal. 
The capacity dimension of the typical snapshot attractor 
of the random Henon case and the random delay case is 
1.14 and 1 respectively, as expected from the observation 
that the former attractor is fractal while the latter is a 
one-dimensional curve. 

Moreover, we calculated also the information dimen- 
sion Di of the four attractors shown in Figdl Results 
are given in Tab.l. Notice that the information dimen- 
sion of the long trajectory attractor for the random delay 
case is close to 1. 

To characterize the dynamical behavior of the system 
we calculate the Lyapunov spectrum of the random delay 
case via the standard method [10]. In general the system 
has one Lyapunov exponent of finite value and the other 
of value — oo due to the dimension collapse. The Kaplan- 
Yorke dimension Dx therefore has the value 1, consistent 
with the obtained value Di 1 for the long trajectory 
attractor as depicted in Tab.l. Furthermore, as shown 
in Figl3] the leading Lyapunov exponent of the system 
can be positive or negative depending on the probabil- 
ity p. The capacity dimension I^o, however, stays con- 
stant in the regime p > 0.2 irrespective of the value of p. 
The smaller value of I^o obtained in the regime p < 0.2 
is a numerical artifact due to an insufficient number of 
phase points used. Simulations show indeed that the ob- 
tained value of 1^0 for p < 02 increases with the number 
of phase points used and the regime of constant I^o ex- 
pands correspondingly. An additional interesting point 
in Fig. [3] is that a randomly temporal switching between 
two nonchaotic dynamics with trivial attractors results in 
a chaotic dynamics holding a fractal attractor (e.g. for 
a = 0.9, k = 0.35 and p = 0.5). Without considering the 
positive Lyapunov exponent corresponding to the exter- 
nal noise ct the fractal attractor of the random delay case 
with a negative Lyapunov exponent may be viewed as a 
strange nonchaotic attractor [14]. 

Dimension collapse. A characteristic feature of the 
random delay system Eq.([T]) is that the length of the de- 
lay time ffuctuates between and 1 randomly. It is there- 
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FIG. 1: Long- trajectory and snapshot attractors, respectively, 
for the random Henon map with b — 0.3, k — 0.05 and p — 0.5 
(a)-(b) and for the random delay system with 6 = 0, /c = 0.3 
and p = 0.5 (c)-(d). In (c) two runs with different initial 
conditions and different realizations of the noise are shown. 
In (d) two attractors at different instants of time are shown. 
The corresponding capacity and information dimensions are 
listed in Tab. I. 
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FIG. 2: Capacity dimension of the snapshot and long- 
trajectory attractors for the random delay map and the ran- 
dom Henon map as presented in[T] 



TABLE I: Calculating the capacity and information dimen- 
sion of the four attractors in Fig. 1 

Do Di 

Fig.l-a 1.95 1.85 
Fig.l-b 1.14 1.06 
Fig.l-c 1.27 1.09 
Fig.l-d 1.00 0.84 



fore convenient for our analysis to partition the evolution 
of the system dynamics into segments having either pure 
delay time or 1, wherein the system dynamics simply 
corresponds to the logistic or the Henon map, respec- 
tively. For simplicity, we will denote them as 0- or 1- seg- 
ment respectively. Whenever the system switches from a 
1-segment to a 0-segment, the phase point representing 
the dynamic evolution of the system will fall immedi- 
ately on the parabola z^+i = 1 — azn and the evolution 
of the system dynamics sticks in such a one-dimensional 
subspace during the 0-segment. This is what we called 
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FIG. 3: (Color online) Lyapunov exponent A (curves) and 
the capacity dimension Do (symbols) of the long- trajectory 
attractor vs. probability p for two random-delay cases with 
5 = 0. Notice the existence of transitions between stable and 
chaotic states. 



dimension collapse. At the end of the 0-segment, the sys- 
tem enters a 1-segment and the dynamics is now just the 
Henon map without randomness until the appearance of 
the next 0-segment. It is obvious that the target of the 
collapses is not influenced by the history of the system 
dynamics or the delay variations. The dynamic evolu- 
tion of the random delay system Eq.([T]) can therefore be 
viewed as the re-injections to the parabola during the 

0- segments and the iterations of the parabola under the 
Henon dynamics during the 1-segments. The long trajec- 
tory attractor of the system can thus be roughly viewed 
as the union of the parabola and its n-fold image under 
the Henon dynamics where n is the possible length of the 

1- segments and n = 1, 2, • • • , +oo. 

An expression of the natural measure p* of the attrac- 
tor can be worked out accordingly. Based on the above 
discussions one may decompose p* in the following way 



Plo 



n=0 



oo 

E 

n=l 



PLn 



(2) 



where pLo is the part of the natural measure living on the 
parabola given by the logistic map and pLn represents the 
part of the nature measure whose support is the n-fold 
iteration of the parabola under the Henon dynamics. On 
the other hand a self-consistent equation for p* can be 
written down as 



p* =pHp* + il-p)Lp* 



(3) 



which is a Frobenius-Perron like equation for our random 
system Eq.([T]). Here H and L denote the evolution oper- 
ators for the density of the logistic and the Henon map 
respectively. Inserting (j2j) in (j3j) reads 

oo oo 
PLO + ^ PLn = (1 - P)LP" + P ^ HpLn^ (4) 



n=0 



Identifying terms on two sides of @ which have the same 
phase space support leads to 

PLO = (i-p)V; (5) 

PLn = P^H'^Pm for n > 0. (6) 
Inserting (|6]) in ([2]) generates a new expression of p* 

oo 

p* = Y,p^H^PLo. (7) 

n=0 

As a consistence check one can insert ([5]) in ([7j), it reads 

CO 

p* = {1-p)J2p"H"U* (8) 
= {l-p){l-pH)-'Lp*, (9) 



which is just a simple reformulation of the Frobenius- 
Perron equation given in (j3|). 
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Eq.(j5j) indicates that p^o is the projection of the nat- 
ural measure p* on the parabola of the logistic map. It 
therefore has a continuous density on this one dimen- 
sional subject for the used parameter. As discussed the 
support of pLn is the n-fold iteration of the parabola un- 
der the action of the Henon map and it approaches the 
Henon attractor as n goes to +00. Dimensional collapse 
means that the dimension of pLn is smaller than the com- 
plete dimension of the phase space which in consequence 
implies that the iterations of pLn under the action of the 
Henon map are overlapped only at measure zero points. 
One would thus expect from Eq.® that the capacity di- 
mension D{) of the (long trajectory) attractor is smaller 
than the dimension of the phase space while the informa- 
tion dimension Di should be determined by pLO and has 
the value 1. This explains the appearance of a fractal at- 
tractor in our random delay system Eq. ([T]) and the above 
obtained values of Di of the corresponding attractors. It 
is necessary to point out that the capacity dimension I^o 
of the (long trajectory) attractor is in general different 
from that of the corresponding Henon attractor (see the 
case a = 0.9, h = 0.35 in Fig. [3]), although they are very 
close if the corresponding Henon dynamics is chaotic. 

Moreover as can be seen from Eq. (|6]) , while other Renyi 
dimensions may vary with the probability p the capac- 
ity and information dimensions of the long- trajectory at- 
tractor are independent of p (see Fig. [3] for numerical 
confirmation). 

Note that the random occurrence of segments of lo- 
gistic map dynamics in our random delay system inter- 
rupts the previous Henon map evolution and starts a new 
Henon evolution with a different point on the parabola 
of the logistic map as the initial condition. The random- 
ness of the delay variation can only determine when the 
reinjection happens and prepare a new initial condition 
for the next Henon evolution. This is obviously irrele- 



vant for the global geometrical properties of the long- 
trajectory attractor. Therefore different realizations of 
external noise Ct would lead to the same long trajectory 
attractor as also indicated by Eq.© (see Fig.l (c)). The 
dynamical evolution of phase point surely depends on 
the realization of the external noise, see for example the 
variation of the shape of the snapshot attractor. In this 
sense the snapshot and long-trajectory attractors contain 
different information of the system dynamics. This is as- 
sumed due to the nonstationary of the evolution rule of 
random systems. 

As can be seen from above analysis the collapse of 
the phase space dimension works as a new mechanism 
of generating fractal attractor in random delay systems 
like ([1]) . The collapse ensures that all segments of Henon 
evolutions starts synchronously from the same low di- 
mensional subspace, the parabola of the logistic map. 
Such a synchronization prevents the randomness to blur 
the fine structure of the attractors generated from the 
Henon evolutions. 

Summary. As a first step towards a deep understand- 
ing of the behavior of systems with fluctuating delay 
time, we studied a simple map system with either one 
or two step memory of the past states. Numerical sim- 
ulations and supplementary arguments showed that the 
dimension collapse of the random delay system has some 
very interesting consequences in the system dynamics for 
instance the fractal properties of the attractor. 
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